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ABSTRACT 



Context. T Tauri stars have X-ray luminosities ranging from Lx = 10 28 — 10 32 ergs _1 . These luminosities are similar 
to UV luminosities (Luv ~ 10 30 — 10 31 ergs _1 ) and therefore X-rays are expected to affect the physics and chemistry 
of the upper layers of their surrounding protoplanetary disks. 

Aims. The effects and importance of X-rays on the chemical and hydrostatic structure of protoplanetary disks are 
investigated, species tracing X-ray irradiation (for Lx > 10 29 ergs _1 ) are identified and predictions for [Ol], [Cn] and 
[N n] fine structure line fluxes are provided. 

Methods. We have implemented X-ray physics and chemistry into the chemo-physical disk code ProDiMo. We include 
Coulomb heating and H2 ionization as heating processes and primary and secondary ionization due to X-rays in the 
chemistry. 

Results. X-rays heat up the gas causing it to expand in the optically thin surface layers. Neutral molecular species are 
not much affected in their abundance and spatial distribution, but charged species such as N + , OH + , H20 + and HaO + 
show enhanced abundances in the disk surface. 

Conclusions. Coulomb heating by X-rays changes the vertical structure of the disk, yielding temperatures of ~ 8 000 
K out to distances of 50 AU. The chemical structure is altered by the high electron abundance in the gas in the disk 
surface, causing an efficient ion-molecule chemistry. The products of this, OH + , H2 + and H30 + , are of great interest 
for observations of low-mass young stellar objects with the Herschel Space Observatory. [0 1] (at 63 and 145 /im) and 
[Cn] (at 158pm) fine structure emission are only affected for Lx > 10 30 erg s _1 . 

Key words, protoplanetary disks: X-ray impact - disk structure - ionized molecules 



1. Introduction 



Youn g stellar objects emit X-ray radiation (|Kovama et all 
I1994D. which is associated with magnetic coronal processes 
dFlaccomio et all 120051: IWolk et all 120051: llmanishi etafl 



1200 ID iets (IPravdo et al 
rays (jGiidel et all 
gion between the star 



200l¥ mostly producing soft X- 



or outflows launched in the re- 
and the disk as predicted in the 



X-wind model ()Shu et all I2000D . Non- accreting weak-line 
T Tauri stars (WTTS) are found to be more X-ray lu- 
minous than classica l accr e ting T Tauri stars (C TTS) 
(jStelzer fc Neuhauserl . 120011: iFlaccomio etHl I2003l 2005). 
because CTTS absorb part o f the X-rays in the accretion 
column (jGregorv et all l2007fl . 

The effects of stellar X-rays on protoplanetary disks 
around T Tauri stars have been studied by many groups: 
iGlassgold et al.l (|2004 120071 ) found X-rays to be important 
in heating the surface layers of the disk and estimated the 
strength of the fine-structure emission lines of [Neil] and 
[Ne ill] , which arise from the warm disk atm osphere ex- 
posed to X-ray radiation. iNomura et al.l (|2007h focused on 
H 2 level populations and line emission, finding X-rays to 
control H 2 level p opulation pump i ng to g ether with UV 
Following work of iMeiierink et all (|2008l ). lErcolano et al.l 



and lGlasseold et all (|2009l ) studied atomic-line diag- 
nostics of the inner regions of protoplanetary disks and for- 
mation of wate r in th e warm disk atmosphere, respectively. 
lErcolano et al.l (|2008l ) and subsequent work focused on the 
X-ray role in c ausing photo-evapor a tion of the surface lay- 
ers of the disk. IGorti fc Hollenbachl (|2004 I2008D include X- 
rays to predict [Aril, [Neil], [Fel], [Si], [Fen] and [Sill] as 
good i ndicators of gas physics in the disk. lWoods fc Willacvl 
(|2009f ) stud ied carbon isotope fr a ctiona tion using the same 
method as IGorti fc Hollenbachl (l2004|) to calculate X-ray 
ionization rates. Recent work of lHenning et al.l (|2010l) stud- 
ied the role of UV and X-rays on C2H column densities and 
excitation conditions in disks around T Tauri stars. 



The scope of this paper is to take a step back and first 
perform a study of the impact of X-rays on the 2D disk 
structure, on molecular ionized species and on observational 
diagnostics such as [0 1] , [C 11] and [N 11] . These molecular 
species were not studied in the above listed previous pa- 
pers, but are of fundamental interest in the l ight of their 
recent detections with the Herschel satellite (jBenz et all 
120101 : iBruderer et all |2010T > . We undertake an exploratory 
study in this letter of the relative effects of UV and X-rays 
on the protoplanetary disks hydrostatic, thermal and chem- 
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Table 1. Examples of molecular dissociation and cross sec- 
tion calculation: CO dissociation can end up in three dif- 
ferent channels which are summed up to give the total CO 
cross section, CH is assumed to follow only one path. 



Table 2. Parameters used in the models: the X-ray lumi- 
nosity values correspond to Models from 1 to 5. 



Molecule 


Product 1 


Product 2 


Cross Section 


CO 







o~co = l/3<7c + 




C+ 


0+ 


1/3(0.5<7 C +0.5ao ) + 




c 


Q2+ 


l/3ao 


CH 


c 2 + 


H 


o"CH = crc 



Molecules included: H 2 , CH, NH, OH, CN, CO, N 2 , SiH, 
NO, 2 , SiO, CH 2 , NH 2 , H 2 0, HCN, C0 2 



ical structure using a series of X-ray models with different 
Lx (X-ray luminosity). 



2. Model 



The 



chemo-physica l disk 
(jWoitke et all 120091: iKamp et all 



mod el ing code 
120101) has 



ProDiMo 
been up- 



dated with X-ray physics and chemistry. 



2.1. Input spectrum 

The incident stellar spectrum used in IWoitke et al.l ([20091 ) 
is composed of a solar model with T fi ff = 5 800 K and 
the c romospheric fluxes of HD 129 333 (jDorren fe Guinanl . 
119941 ). The UV luminosity between 91.2 and 205 nm is 
L\jy ~ 4 • 10 31 erg s . Different models fo r the X-ray in- 
put sp ectrum are adopted in the literature: iNomura et al.l 
(2007) use a two temperature thin therma l plasma model 
to fit the observed TW Hydrae spectrum; lErcolano et al.l 
(|2008|) built synthetic coronal spectra using line and con- 
tinuum emiss ivities from the CH I ANTI compilation of 
atomic data; iGorti fc Hollenbachl (|2008l) use a thermal 



blackbody. We p l ace th e X-ray source on the star and follow 
IClassgold et"aLl (|1997h . Ilgea fc GlassgoTdl (|1999h in using 
an analytic input spectrum to describe the Bremsstrahlung 
emission from an isothermal plasma: F(E) oc 1/E ■ 
exp(— E/kTx) where E is the photon energy between 0.1 
and 100 keV, kTx — 1 keV, and Tx is the plasma tem- 
perature. The X-ray input s pectrum is added to the input 
spectrum shown in Fig. 2 in IWoitke et al.l (|2009f ). 



2.2. X-Ray chemistry 

The following paragraphs describe the implementation of 
primary and secondary ionization in the code. 41 primary 
ionization reactions and 16 secondary ionization reactions 
were added in the chemical network. 



Primary Ionization 

The pr imary rates for X-ray abso rption are calculated fol- 
lowing iMeiierink fc Spaansf (120051) (Appendix D.3. 1 ). Th e 
cross sections are taken from IVerner fc Yakovlevl dl99l . 
Since X-rays are likely absorbed in the K-shells, we assume 
that every X-ray absorption leads to a single ionization for 
H, He, Si, CI and a d o uble i onization for C, N, O, S and Fe 
(jMeiierink fc Spaansl . [2005) . We take into account molec- 
ular X-ray absorption, which leads always to dissociation, 
of the species. Table Q] lists two examples of how the to- 
tal dissociation rates are calculated. When the difference in 



Quantity 


Symbol 


Value 


Stellar mass 


M t 


1 M 


Stellar luminosity 




1 L e 


Disk mass 


M disk 


0.01 M 


X-ray luminosity (0.1-50 keV) 


L x 


0, 10 29 , 10 30 
10 31 , 10 32 


Inner disk radius 


Rin 


0.5 AU 


Outer disk radius 


Rout 


500 AU 


Surface density power law index 


e 


1.5 


Min. dust particle size 




0.1 /an 


Max. dust particle size 




10 /im 


Dust size distribution power index 


f^pow 


2.5 



weight of the elements that form the molecule is large (e.g. 
CH, OH, H2O etc.), we use the cross section of the heavier 
element, when the weight is comparable (e.g. CO, CN, NO 
etc.) we combine the cross sections of the single elements. 

Secondary Ionization 

The fast electrons generated by the X-ray absorption can 
further ionize the gas. The rates are computed as shown 
in Meijerink & Spaa ns (2005, Appendix D.3.2) using ex- 
perimental data from lLennon et all (|l988l ). The volumetric 



rates are a function of the local chemistry, i.e. they depend 
on the species densities jih, "•h 2 ; n <n> an d n c (where uh 
is the hydrogen atom density and n<H> is the hydrogen 
nuclei density). This has to be taken into account when 
solving the chemical equilibrium. The additional entries in 
the chemical Jacobian that come from the X-ray secondary 
ionization reactions are implemented following Sect. 5.6 in 
IWoitke et all (l2009h . 

2.3. X-Ray heating 



We added Coulomb heating (|Dalgarno et all I1999D and H 2 
ionization heating (Meijerink & Spa ans, 2005, Appendi x 
B.l) to the heating proces ses listed in IWoitke et al.l J2009). 



The iMalonev et alT(11996D heating rate only holds for gas 
with low x*. while iDalearno et al.l ([19991 ) present the more 



general case for different high and low ionization gases 
(atomic/molecular). The heating efficiency increases by a 
factor 7-8 in a highly ionized ato mic gas and a factor 2 
in a highly ionized molecular gas (jShull fc van Steenberd . 
1985). 



2.4. Parameter space 

We computed five different models with increasing Lx to 
compare the results with the UV only case (Model 1) for a 
disk surrounding a T Tauri star (Table [5]). In Model 1 
X- rays are swit c hed o ff. This is the model as presented 
111 IWoitke et all (2009). For more details, we refer to the 
original paper. X-ray models (Model 2 — 5) have the same 
parameters as Model 1, but include stellar X-rays (Lx = 
10 29 - lO^ergs" 1 ). We added CI, C1+ and double ionized 
species (C 2+ , N 2 +, 2+ , S 2+ , Fe 2+ ) to the chemical net- 
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1 10 100 1 10 100 1 10 100 

r [AU] r [AU] r [AU] 

Fig. 1. Gas temperature distribution: Model 1 (UV only) in the left panel, X-ray models with Lx = 10 30 erg s _1 (Model 3), 
Lx = 10 32 ergs _1 (Model 5) in the middle and right panel respectively Contour lines are over-plotted for r ra d = 1 (radial 
dust optical depth at 550 nm, dashed line), r vcr = 1 (vertical dust optical depth at 550 nm, dotted line) and X-ray optical 
depth of one at 1 keV (dot-dashed line). The relative position of these two depth depend strongly on the assumed dust 
properties (see Table [5]). The solid line corresponds to T gas = 8 000 K. 




1 10 100 1 10 100 1 10 100 

r [AU] r [AU] r [AU] 

Fig. 2. Total hydrogen nuclei density n<H> distribution: Model 1 (UV only) in the left panel, X-ray models with Lx = 

10 30 ergs _1 (Model 3), Lx — 10 32 ergs^ 1 (Model 5) in the middle and right panel respectively. The dotted contour line 

shows where the ratio between Coulomb heating and photo-electric plus PAH heating equals 1. The solid line circles the 

hot water reservoir (Th 2 o > 200 K). 



work. The X-ray emission is treated as a point source at 
the location of the star. 

3. Results 

We analyze the results obtained from the five models de- 
scribed above. We first present the thermal and hydrostatic 
structure of the various models. Then we describe the im- 
pact of X-rays on the chemistry and fine structure line emis- 
sion, focusing on the main differences between the UV only 
model and the combined UV+X-ray models. 

Temperature and density structure 

The low density (n<n> < 10 s cm -3 ) surface layers di- 
rectly exposed to X-ray radiation are heated efficiently via 
Coulomb heating: in all X-ray models temperatures reach 
up to ^8 000 K (Fig.[T]). The extension of this high temper- 
ature region increases with Lx- it reaches 2 AU in Model 2 
and 100 AU in Model 3 and Model 4. In Model 5 the tem- 



perature is slightly lower for r > 50 AU but the whole disk 
is warmer (Fig. [I]). X-rays dominate the thermal balance 
only in the upper layers down to (X/UV)heat = 1- Beyond 
(X/UV)hoat = 1, UV heating takes over and becomes the 
main heating process. Deeper in the midplane, both X-ray 
and UV heating are negligible, but their ratio is inverted 
again and harder X-rays (E > 10 keV) deposit more energy 
than UV radiation (Fig. EJ. 

The inner wall of the disk is always directly illumi- 
nated by the stellar radiation causing the inner rim (0.5 
AU < r < 0.7 AU) to have a large scale height and hence 
vertical extent. In the UV model the shadowed disk atmo- 
sphere beyond the inner rim has a lower vertical scale height 
since radiation cannot penetrate efficiently and sustain high 
gas temperature. Only radiation that impinges the disk at 
higher angles (z/r > 0.6) and scattered radiation reach the 
outer disk causing the temperature to increase and con- 
sequently the second bump at ~ 9 AU. When X-rays are 
switched on, the temperature increase in the surface layers 
affects directly the vertical density distribution there since 
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Fig. 3. Oxygen line flux emission at 63 and 157 /im, ionized 
carbon line flux emission at 157 /im and ionized nitrogen 
line flux emission at 205 fim as a function of the X-ray 
luminosity. All fluxes are compiled for a distance of 140 pc. 
The correlation between N n and Lx is clear, while O I and 
Cll show no correlation for Lx < 10 30 erg s" 1 . 



Fig. 4. Column density of OH + (upper left panel), H20 + 
(upper right panel), H30 + (lower left panel) and N + (lower 
right panel) at three different radii (3, 10, 100 AU) versus 
X-ray luminosity. In each panel the UV only model is indi- 
cated with the "UV" label on the x-axis. 



we keep the radial surface density profile fixed. The inner 
rim of the disk (r < 1 AU and z/r > 0.1) is more puffed up 
(Fig. [21). The second bump at ~ 9 AU is more vertically ex- 
tended as well; in Model 5 it even "merges" with the inner 
rim and flares as a whole (Fig. [3J third panel). 

Chemistry 

The abundances of neutral molecules such as CO and OH 
and H2O generally change by less than an order of magni- 
tude over the entire modeling space. On the other hand, 
the ion chemistry in the disk surface is strongly changed: 
the electron abundance in the upper layers increases by 
a factor 10 4 . This leads to an enhancement of OH + , 
H2O" 1 ", H3O" 1 " and N + abundances. Secondary ionization 
of H2 increases the Hj abundances. H^ collisions with 
H2 produce Hj. Subsequent collisions of with O, OH 
and H2O lead to the production respectively of OH + , 
H20 + and H3O" 1 ". Fig. 2] shows how the vertical column 
densities of these species change with Lx for three different 
disk radii: 3, 10 and 100 AU. OH+ and H 2 0+ respond 
more efficiently than HaO + to the X-ray radiation as the 
total mass ratio M(OH+)/M(H 2 0+) is only doubled from 
Model 1 to Model 5, while M(OH + )/M(H 3 0+) increases by 
more than a factor 5. 

OH + : OH + formation is enhanced approximately propor- 
tional to the X-ray luminosity. Its total mass increases 
from 8 x 1O~ 15 M in Model 1 by a factor 3.6 in Model 3 
and about 500 times in Model 5 (Table EJ). 

As the X-ray luminosity increases, OH + formation is 
enhanced and pushed outwards, where the disk is colder. 
The column density of OH + increases substantially beyond 
r > 10 AU in all models. It reaches ~ 100 times its initial 
value (Model 1) in Model 3 at 100 AU (Fig.g]). 

H20 + : The H20 + mass increases by 3 times (factor 
of 100) between Model 1 and Model 3 (Model 5). The 



column density again increases toward the outer part of 
the disk. Model 3 has ~ 10 times higher column density 
than Model 1 at 10 and 100 AU. As in the case of OH+, 
H 2 + is pushed toward the outer part of the disk with 
increasing X-ray luminosity. 

HaO + : The H30 + mass is slightly less affected by 
X-rays: it is 2 times higher than in the UV case in model 3 
and 10 times in model 4 . In model 3, model 4 and model 5, 
H 3 + is also clearly pushed to larger radii and lower 
temperatures. 

N + : We find the most extreme mass increase for N + 
in all X-ray models (Table 2]) and at all radii: in model 3 it 
is 1 000 times higher than in the UV only case. 

In models 4 and 5, the N + mass increases over- 
proportionally in the outer part of the disk (Fig. |4j right 
lower panel) . Figure [5] shows the column density ratio of 
the X-ray dominated N + to the UV dominated C + . With 
increasing Lx, the N + column density is enhanced by two 
(outer disk) to six (inner disk) orders of magnitude. 

Line Predictions 

In table|3jwe list the flux prediction of [0 1] 63 /im, [0 1] 145 
/im, [C 11] 157 /im and [N 11] 205 /im for all the models. Fig. 
[3] shows their values as a function of the X-ray luminos- 
ity. The [N 11] line flux shows the strongest correlation with 
the X-ray flux, increasing by 5 orders of magnitude from 
model 1 to model 3. The oxygen fine structure line fluxes 
and the [C 11] line flux are constant from model 1 to 3 (fig. 
[3]), due to the fact that they originate in the UV heated 
layer. Only in model 4 and 5 the line fluxes increase with 
L x . 
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Table 3. Predicted line fluxes of [Oi] at 63 /im, [Cn] at 157 /im and [Nil] at 205 /im for all the models expressed in 
W/m 2 . All fluxes are computed for a distance of 140 pc. 



Line 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


[0 i] 63 fim 
[0 I] 145 /urn 
[C n] 157 fim 
[N ii] 205 fim 


4.8(-17) 
1.5(-18) 
1.5(-17) 
3.6(-26) 


4.6(-17) 
1.3(-18) 
1.5(-17) 
5.2(-23) 


4.8(-17) 
1.4(-18) 
1.7(-17) 
4.6(-22) 


8.3(-17) 
2.3(-18) 
3.0(-17) 
4.1(-21) 


3.6(-16) 
1.2(-17) 
6.1(-17) 
l.l(-20) 



Table 4. Species masses given in units of solar mass in all the models 



Species 


UV Only 


L = lO^ergs -1 


L = 10 au ergs -1 


L = lO^ergs- 1 


L = 10^ ergs" 1 


H 


5.0(-05) 


5.0(-05) 


5.2(-05) 


7.1(-05) 


1.4(-04) 


H 2 


7.6(-03) 


7.6(-03) 


7.6(-03) 


7.6(-03) 


7.5(-03) 


C 


3.0(-07) 


3.0(-07) 


3.1(-07) 


3.4(-07) 


4.8(-07) 


c+ 


1.3(-07) 


1.3(-07) 


1.3(-07) 


1.4(-07) 


1.8(-07) 


o 


1.4(-05) 


1.4(-05) 


1.4(-05) 


1.4(-05) 


1.4(-05) 


CO 


9.5(-06) 


9.4(-06) 


9.3(-06) 


9.0(-06) 


9.0(-06) 


OH 


3.4(-ll) 


3.4(-ll) 


3.7(-ll) 


7.0(-ll) 


2.8(-10) 


H 2 


7.7(-07) 


4.7(-08) 


3.6(-08) 


2.7(-08) 


3.3(-08) 


OH+ 


7.9(-15) 


1.1(— 14) 


2.9(-14) 


4.1(-13) 


4.6(-12) 


H 2 0+ 


1.9(— 15) 


2.1(-15) 


5.3(-15) 


6.2(-14) 


4.2(-13) 


H 3 0+ 


3.4(-15) 


3.7(-15) 


7.3(-15) 


4.8(-14) 


3.5(-13) 


N+ 


2.9(-15) 


8.7(-13) 


3.9(-12) 


9.9(-12) 


3.6(-ll) 



4. Discussions and conclusions 

We find that X-rays mainly impact the surface layers of pro- 
toplanetary disks. The hydrostatic structure of these layers 
is significantly changed: X-ray radiation is mostly absorbed 
in the tenuous layers (z/r > 0.2 cm -3 ) increasing the tem- 
perature up to ~ 8 000 K; at those temperatures cooling by 
[O I], Fe II and Lya balance the X-ray heating. This causes 
the density structure to flare stronger than in the UV only 
case (Fig. (2). 

The general resul ts from iNomura et al.l (|2007l ) and 
iGlasseold et al.l (|2004l ) are qualitatively confirmed, show- 
ing that the inner rim and the surface layers are dom- 
inated by X-rays. In addition we find that more distant 
parts, up to 100 AU, are affected as well (Fig. 2]). We find 
that the size of the high temperature X-ray driven region 
increases with Lx- This is due both to the fact that our 
model uses the iDalearno et al.l (|l999t ) prescription for the 
X-ray heating and to o ur choice of the input spectrum. The 
IDalearno et al.l (|1999f) heating efficiencies lead to ~ 7 — 8 
times higher rates in region s with high electron fraction, 
compared to iMalonev et al.l (|1996[ ). Furthermore our in- 
put spectrum has higher photon flux in the range .1-0.3 
keV compared with e.g.. iGorti fe H ollcnbach (2008) and 
INomura et al.l (|2007l ). These soft X-ray photons are ab- 
sorbed at low vertical column densities causing high e~ den- 
sities there and thus are more efficie nt in heating the u pper 
layers. This effect was also noted in lGorti et al.l ((2009). 

IGorti fe Hollenbachl (|2008t) found X-rays to be the dom- 
inant heating process in slightly deeper layers than we do. 
This is most likely due to the fact that we include UV 
scattering from dust, which leads to a larger vertica l pene- 
tration of FUV radiation in the disk. INomura et al.1 (|2007h 
also show that FUV penetrates deeper than X-rays, be- 
cause FUV scattering is more efficient than X-ray Compton 
scattering at 1 keV , whic h is not included in our model. 
IGorti fe Hollenbachl (|2008f ) also consider a less steep surface 
density distribution which yields less mass at 1 AU, where 
most of the X-ray radiation is absorbed in our model. 



The X-ray energy deposition does not affect neutral oxy- 
gen fine structure emission at 63 and 145 ^m (Table [3]), 
unless the X-ray luminosity is pushed to Lx > 10 30 ergs -1 , 
where X-ray heating dominates in the regions where the 
emission takes place. Ionized carbon emission at 157 /im 
is only slightly affected by the X-rays in all the mod- 
els, as the line emission is dominated by materi al at R 
> 200 AU (|Woitke et all . l2009t iKamp et all I2010H . where 
the X-ray contribution to the chemistry and thermal bal- 
ance is negligib l e. Ou r val ues are in good agreem ent with 
iMeiierink etall (|2008l ) and lErcolano et al.l (|2008l ) for L x ~ 
10 30 erg s — 1 . However the correlation of th e oxyg en line 
fluxes with Lx, found in the IMeiierink et all (|2008t ) model 
is due to the lack of UV radiation in their model. In our 
model at low Lx, the UV heating provides the gas temper- 
ature to sustain a constant level of [O i] emission. 

We find H2, CO and OH properties in the X-ray mod- 
els largely unvaried with respect to Model 1. The largest 
change occurs for water which mass diminishes by a fac- 
tor 20 in the highest Lx model compared with the UV 
only model. The X-rays enhance the electron abundance in 
the surface layers, where we find x e ~ 0.1 for z/r ~ 0.9 
out to 100 AU. Ionized species abundances are enhanced, 
such as OH+, H 2 0+, H 3 0+ and N+. The molecular ion- 
ized species formation is pushed outwards as Lx increases. 
High abundances of OH + and H20 + are a characteristic 
signature of an X-ray domina t ed re gion as recently dis- 
cussed by Ivan der Werf et aT] (|2010D after the detection 
with Herschel/SPIRE in the ultra- luminous galaxy Mrk 
231. Furthermore this is extremely interesting in the con- 
text of the recent detection of thes e species in a massive 
YSO by Herschel (jBenz et all 12010). The column densities 
that we find from our X-ray irradiated d isks are of the sam e 
order of magnitude as those derived in lBenz et al.l ((2010). 
This suggests that a disk could in principle bring a sub- 
stantial contribution to the line fluxes emitted from these 
molecules, making our result extremely relevant for future 
Herschel studies of low mass YSOs. 
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The abundance of these ions increases especially at 
larger radii (r > 10 AU), while the N + abundance rises 
considerably at all radii. The latter is not observed in the 
UV model due to the high ionization potential of N(I.P.) = 
14.5 eV. We find that the N+ flux emitted at 205 ^m in 
model 3 is ~ 5 x 10~ 22 W/m 2 , coming for 90% from the 
upper layers beyond 10 AU. This is much smaller than 
the current Herschel sensitivity limit of 5 x 10~ 18 W/m 2 
(HIFI, 5 a in 1 hour) and even below current sensitivity es- 
timates for SPICA/SAFARI (10~ 19 W/m 2 , 5 a in 1 hour). 
We only consider X-ray ionization, but the flux is unlikely 
to change if EUV is present. Additional ionization occurs 
only at smaller rad ii, because of the low pene tration depth 
of EUV radiation (jHollenbach fc Gortil . [20091) . The molec- 
ular ionized species formation is pushed toward the outer 
disk as Lx increases. 



5. Outlook 

The predicted N + line fluxes are too low to be observed with 
current instruments. In contrast, the [O i] 63 /xm and [C n] 
157 /im lines are observable in the presence of UV luminosi- 
ties of ~ 10 31 erg s _1 and/or strong X-rays (Lx > 10 30 erg 
s -1 ). In the context of the current Herschel observatory, the 
OH + , H20 + and H3CV lines are extremely interesting as 
they are likely observable with the HIFI, PACS and SPIRE 
instruments. These hydride ions and the already observed 
[Ne II] and [Ar II] line fluxes, will be addressed in the sub- 
sequent paper (Aresu et al., in preparation). We will also 
investigate a larger parameter space and the effect of the 
change in disk structure - especially the higher inner rim 
and increased flaring - on the overall SED. Future work also 
includes the creatio n of a grid of X-ray models following the 
IWoitke et afl <|2010h approach. 
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Fig. 5. Column density ratio between the X-ray dominated 
N + and the UV dominated C + . The latter does not change 
in the X-ray models with respect to the UV model 
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